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We propose a method for measuring spin Hamiltonians and apply it to the spin- 1/2 Heisenberg 
antiferromagnet CS2CUCI4, which shows a 2D fractionalized RVB state at low fields. By applying 
strong fields we fully align the spin moment of CS2CUCI4 transforming it into an effective ferromag- 
net. In this phase the excitations are conventional magnons and their dispersion relation measured 
using neutron scattering give the exchange couplings directly, which are found to form an anisotropic 
triangular lattice with small Dzyaloshinskii-Moriya terms. Using the field to control the excitations 
we observe Rose condensation of magnons into an ordered ground state. 

PACS numbers: 75.10.Jm, 75.45.+j, 75.40.Gb, 05.30.Pr 



Understanding strongly correlated physics poses 
formidable mathematical difficulties and in only a few 
exceptional cases has the full many-body quantum prob- 
lem been solved. Although theory takes the Hamiltonian 
(7i) as its starting point linking experimental data to this 
is often not possible. A method for measuring TL directly 
would bridge this gap to theoretical approaches and in 
addition reveal the essential ingredients from which ex- 
otic quantum states emerge. Motivated by this wc com- 
bine neutron scattering with high magnetic fields and 
make just such a determination of Ti. taking the remark- 
able quantum magnet CS2CUCI4 as a subject. We base 
our approach on overcoming spin couplings using large 
fields thus transforming the system into an effective fer- 
romagnet, an easily solvable state. In addition we ex- 
plore how the ordered ground state evolves with lower- 
ing field and interpret the results in the framework of 
Bose-Einstein condensation (BEC) of magnons. 

The insulating magnet CS2CUCI4 is an ideal subject 
for two reasons: First, its relatively weak (~ 4 K) cou- 
plings can be overcome by current fields (at 8.44 T), and 
second, it shows highly unusual strongly correlated prop- 
erties [0 . Among the most fascinating are a low-field dy- 
namics dominated by 2D highly dispersive continua char- 
acteristic of fractionalization of spinwaves into spin- 1/2 
spinons, exceptionally strong quantum renormalizations, 
and an unexplained T — Q disordered phase induced by 
weak fields along h and c. Although Anderson first pro- 
posed a 2D fractionalized state in 1973 (the resonating 
valence bond state), the essential conditions for its exis- 
tance have remained highly contentious In light of 
this establishing what the special ingredients are in the 
Hamiltonian of CS2CUCI4 is therefore very important. 

The origins of strongly correlated phases lie in the un- 
certainty principle. For quantum magnets uncertainty 



is embedded in the noncommutation of the spin vector 
components S = {S^ , 3^,3^} and the "true" direction of 
S cannot be known. For spins Sn on a lattice R coupled 
by the Heisenberg exchange Hamiltonian 

^ = ^ JsSr ■ Sr+s - guBBS^ (1) 

the energy depends simultaneously on all three non- 
commuting components of each Sn {S is a vector be- 
tween sites and the last term an attendant magnetic 
field). Quantum uncertainty appears as a kinetic term 
i'^R'^R+s + ^R'^R+s) action on the spins which is 

most extreme for spin-1/2 where it fiips pairs of spins e.g. 
Ti to It, and the magnet fluctuates between many spin 
configurations. Semiclassically this kinetic action corre- 
lates particle motions (and creation) with others and can 
be so strong that new phases emerge as in CS2CUCI4. 

When large enough, the field B in (1) prevails over 
the exchanges and the unique situation arises where the 
ground state of Ti is known and the one-particle ex- 
cited states are exactly solved. The ground state con- 
sists of all spins up, which we denote V&o = |0), Eq = 
-NgfisBS -\- N/2J2sJs3'^, which is indeed that of a 
ferromagnet. There are N orthonormal states with a sin- 
gle spin fiip Vk = S^/V^\0) corresponding to all sites 
R. When Ti acts on tpR, it generates only other such 
one-spin-flip states because the total spin S'f. = J2r '^r 
is a constant-of-the-motion for H. Because the Hamil- 
tonian is invariant upon translation plane-wave states 
are diagonal: ipk = N~^^^^jj^exp{—ik ■ R)ipR where 
TLipk = E{k)ipi^. The kinetic term causes hopping of 
these spin-flips through the lattice and the energy eigen- 
values (for spins 5*= 1/2) are 

Eik) =Eo+ g^iBB - Jo + Jk. Jk = \Y. ■^^'"''''^ (2) 
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so that the one-spin-flip excitations disperse relative to 
the ground state with the relation, Tiujk = E(k) — Eq = 
gUsB — Jq + Jky which is a constant term plus J^, 
the Fourier transformed exchange couplings. These ex- 
citations are the familiar quantized harmonic spin-wave 
modes, magnons, which carry AS^ = — 1 and have Bose 
statistics. Since neutrons are spin-1/2 particles they scat- 
ter only by changing the total spin by AS' = 0, ±1. 
For a system prepared in the fully-aligned state |0) neu- 
trons can scatter inelastically only by exciting a single 
magnon through the matrix element KV'fcl'S'fe |0)P where 
S7 = A^^^/^ ^^exp(— ife • R)S^. Jk and so Jg of Eq. 
(|l|) can then be found from the measured dispersion huk- 
The crystal structure of CS2CUCI4 is orthorhombic 
(Pnma) with lattice parameters at 0.3 K of a=9.65 A, 
6=7.48 A, and c=12.35 A. The magnetic S = 1/2 Cu2+ 
ions are situated within distorted CuCl4~ tetrahedra. 
Layers {be plane) of these tetrahedra are separated by 
Cs"*" ions and are stacked with an offset giving the struc- 
ture illustrated in Fig. |^(a). The strongly correlated 
physics derives from the "isosceles" triangular lattice ar- 
rangement of spins in the layers with antiferromagnetic 
exchange paths J = J[o,i.o] ^^'^ J' = "^[o, 1/2, ±1/2] • The 
triangular geometry allows a large configuration space 
for fluctuations and is presumably crucial to fractional- 
ization. 

To measure hujk the V2 cold-neutron triple-axis spec- 
trometer at the BER-II reactor at HMI in Berlin was 
used. A large (3.6 g) high-quality single crystal of 
CS2CUCI4 was mounted in the (0,fc,Z) scattering plane 
on a dilution refrigerator insert with base temperature 
of 50 mK. The VM-1 cryomagnet provided fields up to 
14.5 T along a. The spectrometer was configured with 
a vertically-focused monochromator (PG002) and a hori- 
zontally focused PG002 analyzer to select scattered neu- 
trons with fixed fc/ = 1.2 or 1.35 A~^. 

A magnetic field of 12 T, much larger than the sat- 
uration field (8.44 T), was used to open a significant 
energy gap of 0.435(8) meV to the first excited states. 
Temperatures below 200 mK ensured that the thermally 
introduced population of spin flips was less than 1 per 
10^^ spins. No magnetostructural distortions were ob- 
served and the origin of superexchange in high-energy 
electronic bonds means that the coupling constants are 
unperturbed by the field. Only one magnon scattering 
events were observed and their energy and wavevector 
dependence mapped out. Fig. ||(c) shows a typical scan. 
Two resolution limited peaks are seen separated by a 
small energy of 0.084(2) meV; this splitting is due to an 
additional anisotropy as explained below. 

The measured one-magnon dispersion relations are 
graphed in Fig. ^(a). The considerable dispersion along 
both [OfcO] and [OOZ] in the be plane indicates strong 
2D character. The overall dispersion follows hujk with 
Jk = J cos{2TTk) + 2J' coa^nk) cos(7r/) where J = 0.374(5) 
meV and J' = 0.128(5) meV, the couplings in Fig. |^(a) 



TABLE I: Hamiltonian parameters {B > Be) (see textWersus 
the quantum renormalized {B = 0) parameters from 



Parameter 


B> Be 


B = 


Renormalization 


J (meV) 


0.374(5) 


0.62(1) 


1.65(5) 


J' (meV) 


0.128(5) 


0.117(9) 


0.91(9) 


J" (meV) 


0.017(2) 






Da (meV) 


0.020(2) 






e (rlu) 


0.053(1) 


0.030(2) 


0.56(2) 



(fe — {h, k, I) is expressed in units of (27r/a, 27r/&, 27r/c)). 
The small splitting into two magnon branches is charac- 
terized well by modified dispersions huj^ = g^sB — Jo + 
Jk ± Dk,Dk = 2DaSin(7rA:)cos(7r/) with Da = 0.020(2) 
meV. This is surprising because whereas Jk is a sum of 
cosine terms, Dk is sinusoidal. The physical meaning of 
this is that a left moving magnon (of a certain type) has 
different energy from a right moving one hio^ ^ ^^-k'y 
the two magnon branches actually cross over at A: = 0, 
hujf — huj'^f,, and such a situation can come about only 
if an exchange with a sense of direction is present. 

Dzyaloshinskii and Moriya (DM) proposed just such 
an exchange interaction many years ago. They showed 
that spin-orbit couplings in the superexchange can gen- 
erate a coupling of the form Dij ■ {Si x Sj). Their in- 
teraction is of higher-order and therefore much weaker 
than Heisenberg exchange and can occur only when the 
superexchange pathways do not have centers of inversion 
which is indeed the case in CS2CUCI4. 

In the ordered structure (below Tn—620 mK and 
-0 = 0) the moments lie almost within the be plane 
which would happen if the most important Dij vec- 
tor in CS2CUCI4 was directed along the a-axis, and the 
observed wavevector dependence, sin(7rA;) cos(7rZ), of Dk 
suggests that the important DM interaction is along the 
same zig-zag bonds as J' in the 2D planes. Considering 
these bonds only, and making the approximation that 
= {±Da, 0, 0) = we obtain using symmetry 

R 

(3) 

where the labels <5i_4 refer to Fig. |l](a) and the ± has 
been introduced because there are two distinct layers 
shown in Fig. ^(b) which are inverted versions of each 
other with DM vectors pointing in opposite directions. 
Like the Heisenberg coupling this DM interaction also 
conserves and plane-wave solutions remain diagonal; 
'H^j^.j'ipk = ^Dkipk where Dk = 2DaSm{Trk) cos{Trl) as 
observed. The DM interaction then explains the observed 
sinusoidal components of hco^ and the fact that there are 
two modes - one for each type of layer. 

The fact that CS2CUCI4 orders three dimensionally 
means that there must be an interaction J^' between lay- 
ers. We introduce operators and 6^ that create the 
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two types of magnons on the different layers. The full 
Hamiltonian with DM and interlayer couplings is: 



hQk — Jq + Dk 



I" 



J" 
■^k 



Klk — Jq — Dk 



a-k 
bk 



where hilk = QIJ-bB — Jq + Jk is the magnon dispersion 

= 0. D 
the new dispersion relations 



for Dk = J'k = 0. Diagonalizing this Hamiltonian gives 



huji. 



T" + 



(4) 



and for the case of interlayer nearest neighbor coupling 
[see Fig. ||(b)] {J'k^J" cos(7r/i)e"'2''"^) the relative inten- 
sity of the two modes is 



1 + cos(27r/C) sin(26'fc) 
1 - cos(27r/C) sin(26'fe)' 



(5) 



with Ok = tan-i[J"cos(7r/t)/(/DjT]J^ + Dfc)] and 
where the total inelastic intensity + is indepen- 
dent of wavevector. Here C=0.34 is the relative off- 
set along c between adjacent layers. Fitting the above 
model (Eqs.(||) and to the data yields the excellent 
fits shown in Figs. p|a)-(d) with the fitted parameters 
listed in the first column of Table I and 5a=2.19(l). The 
total inelastic intensity shown in Fig. ^(b) is nearly in- 
dependent of k as predicted. The relative intensity of 
the two modes (where they could be resolved) is shown 
in Fig. ||(d). We conclude that all other couplings in 
CS2CUCI4 are much smaller. Dipolar energies and g- 
tensor anisotropies are small and neglected here. 

Upon decreasing field the magnon energies reduce by 
the additive Zeeman term gjisB [see Fig. ||(a)]. At the 
critical field i?c=8.44(l) T the gap closes at the dis- 
persion minima r ± Q, Q={0.5+e)b*, e=0.053(l). At 
those wavevectors Bragg peaks appear below Be indicat- 
ing transverse (off-diagonal) long-range order. This order 
is an example of BEC in a dilute gas of magnons induced 
by changing the "chemical potential" \Bc ^ B\ [Q. The 
measured spin order forms an elliptical cone around the 
field direction {Sr) = ±bSb cos Q R + cSc sin Q R + aSa 
(odd/even ± layers contrarotate) where Sb > Sc as il- 
lustrated in Fig. ^(e). In fact this order corresponds ex- 
actly to the simultaneous condensation of contrarotating 
magnons wZq O'f^d oj^q [see Fig. «(d)] with gap closure at 
Be; a mean- field calculation ||^, ^ of this state gives an 
elliptical cone with asymmetry (cos 6'q -|-sin 9q)/ (cos 9q— 
sin (?q)==1. 52(6) in agreement with the observed ratio 
'S'6/S'c=l. 55(10) just below Be- The asymmetry is a com- 
bined effect of interlayer coupling J" and alternation of 
Z)^ between layers and rapidly decreases as the field is 
lowered due to increased inter-particle interactions and 
fluctuations, S'b/5'c=l.l(l) below 7 T. 

The effect of fluctuations and interactions on the order 
as fleld decreases is quantified in Fig. ^(b-c): Fig. ||(b) 
shows the off-diagonal order parameter Sc- Close to Be 



it is described by a power law (solid line) Sc^\Bc — B\^ 
with /3=0.33(3), significantly below the value /3=0.5 ex- 
pected for mean-field (3D) BEC Q]. The magnetization 
Sa obtained from susceptibility measurements 0| is plot- 
ted in the inset of Fig. ||(c). It shows that the boson 
density is not linear versus \Bc — B\ but rather shows 
a deviation that may be logarithmic ||]; and finally the 
wavevector of the condensate Q = 0.5 -|- e{B) is plot- 
ted in Fig. ^(c). Q varies strongly with field indicating 
that magnon-magnon interactions are important even at 
low density and renormalize the condensate wavevector. 
The above features deviate significantly from mean-field 
(3D) behavior Q and could be associated with the 2D 
nature of the magnons. In two dimensions interactions 
can qualitatively change the scaling behavior such as by 
introducing non-linear, log corrections to the magnetiza- 
tion curve lH . 

In summary, we have determined the Hamiltonian of 
the quasi-2D quantum magnet CS2CUCI4 using a new ex- 
perimental method and show that it is a 2D anisotropic 
triangular system. We also measured transverse (off- 
diagonal) order with field below saturation, an example 
of Bose-Einstein condensation of magnons. Our meth- 
ods are general and could be used to reveal exchanges 
and quantum renormalizations for systems as diverse as 
random magnets, quantum antiferromagnets and spin 
glasses. 
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k {iTilb) 

FIG. 1: (a) Magnetic couplings in a 2D triangular layer in 
CS2CUCI4: strong bonds J (heavy lines || h), smaller frustrat- 
ing zig-zag bonds J' (thin lines). Da are the Dzyaloshinskii- 
Moriya (DM) couplings along the zig-zag bonds; the signs 
(g), refer to interactions originating at the central spin Sr, 
see Eq. (|^) . (b) Odd (black) and even (grey) layers are stacked 
successively along a (inter-layer spacing a/2) with an offset 
("=0.34 along c. J" (dashed arrow) is the nearest-neighbor 
inter-layer exchange, (c) (fe*,c*) reciprocal plane showing the 
near-hexagonal Brillouin zones (thin lines) of the triangular 
lattice in (a). Black points are zone centers (t) and grey 
points at r ± <5 mark dispersion gap minima in the saturated 
phase {B > Be) and where Bragg peaks appear below Be- 
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FIG. 2: (a) Magnetic dispersion relations in the saturated 
phase (5=12 T||a, T <0.2 K) along symmetry directions in 
the 2D plane (heavy dashed lines in Fig. |l|(c)). Solid and 
dashed lines are fits to Eq. with parameters in Table 1. (b) 
Observed integrated inelastic intensity compared with predic- 
tions for the fully-polarized eigenstate (solid line), (c) Exci- 
tations lineshape observed along a constant-wavevector scan 
at the minimum gap fc=(0, 1.447,0). Solid line is a fit to Eqs. 
(§^1) convolved with the instrumental resolution (horizontal 
grey bar indicates the full-width-half-maximum of the energy 
resolution) . (d) Relative intensity of the two magnon modes 
compared with Eq. (0) (solid line). 
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FIG. 3: (a) Magnon energies vs. field in the saturated 
phase. Solid lines are fits to a linear behavior as expected for 
AS^ = -1 eigenstates with g-factor pa=2.18(l). (0,0.447,1)=*= 
label the two magnon modes resolved at the minimum gap in 
scans such as in Fig. |2|(c). (b) Amplitude of perpendicular 
ordered moment Sc in the cone phase vs. field. Solid line 
is a power-law fit. (c) Incommensuration {e=Q-0.5) vs. field 
fsolid line is guide to the eye). Inset: magnetization vs. field 
M (r=30 mK) compared with a linear behavior (solid line), 
(d) Superposition of contrarotating magnons tt'^Q |^] of dif- 
ferent amplitudes (large and small circle) gives the elliptical 
order in be plane shown schematically for odd layers in (e) 
(arrows are ordered spins). Even layers have opposite sense 
of rotation. 



